Finite Element Solvers for Biot’s Poroelasticity Equations in Porous Media     

Kadeethum  T.  Lee  S.  Nick  H. M. 《Mathematical Geosciences》2020,52(8):977-1015

We study and compare five different combinations of finite element spaces for approximating the coupled flow and solid deformation system, so-called Biot’s equations. The permeability and porosity fields are heterogeneous and depend on solid displacement and fluid pressure. We provide detailed comparisons among the continuous Galerkin, discontinuous Galerkin, enriched Galerkin, and two types of mixed finite element methods. Several advantages and disadvantages for each of the above techniques are investigated by comparing local mass conservation properties, the accuracy of the flux approximation, number of degrees of freedom (DOF), and wall and CPU times. Three-field formulation methods with fluid velocity as an additional primary variable generally require a larger number of DOF, longer wall and CPU times, and a greater number of iterations in the linear solver in order to converge. The two-field formulation, a combination of continuous and enriched Galerkin function space, requires the fewest DOF among the methods that conserve local mass. Moreover, our results illustrate that three out of the five methods conserve local mass and produce similar flux approximations when conductivity alteration is included. These comparisons of the key performance indicators of different combinations of finite element methods can be utilized to choose the preferred method based on the required accuracy and the available computational resources.

  相似文献   

Wetlands and low-gradient topography are associated with longer hydrologic transit times in Precambrian Shield headwater catchments     

Daniel Lane  Colin P.R. McCarter  Murray Richardson  Chris McConnell  Tim Field  Huaxia Yao  George Arhonditsis  Carl P.J. Mitchell 《水文研究》2020,34(3):598-614

The estimation of hydrologic transit times in a catchment provides insights into the integrated effects of water storage, mixing dynamics, and runoff generation processes. There has been limited effort to estimate transit times in southern boreal Precambrian Shield landscapes, which are characteristically heterogeneous with surface cover including till, thin soils, bedrock outcrops, and depressional wetland features that play contrasting hydrologic roles. This study presents approximately 3.5 years of precipitation and streamflow water isotope data and estimates mean transit times (MTTs) and the young water fraction (p_y) across six small catchments in the Muskoka-Haliburton region of south-central Ontario. The main objectives were to define a typical range of MTTs for headwater catchments in this region and to identify landscape variables that best explain differences in MTTs/p_y using airborne light detection and ranging and digital terrain analysis. Of the transit time distributions, the two parallel linear reservoir and gamma distributions best describe the hydrology of these catchments, particularly because of their ability to capture more extreme changes related to events such as snowmelt. The estimated MTTs, regardless of the modelling approach or distribution used, are positively associated with the percent wetland area and negatively with mean slope in the catchments. In this landscape, low-gradient features such as wetlands increase catchment scale water storage when antecedent conditions are dryer and decrease transit times when there is a moisture surplus, which plausibly explains the increases in MTTs and mean annual runoff from catchments with significant coverage of these landscape features.  相似文献   

Changes in hillslope hydrology in a perched,shallow soil system due to clearcutting and residual biomass removal     

Colin P. R. McCarter  Stephen D. Sebestyen  Susan L. Eggert  Randall K. Kolka  Carl P. J. Mitchell 《水文研究》2020,34(26):5354-5369

Sustainable fuels legislation and volatility in energy prices have put additional pressures on the forestry sector to intensify the harvesting of biomass for “advanced biofuel” production. To better understand how residual biomass removal after harvest affects forest hydrology in relatively low slope terrain, a Before-After-Control-Impact (BACI) study was conducted in the USDA Forest Service's Marcell Experimental Forest, Minnesota, USA. Hydrological measurements were made from 2010–2013 on a forested hillslope that was divided into three treatment blocks, where one block was harvested and residual biomass removed (Biomass Removed), the second was harvested and residual biomass left (Biomass Left), and the last block was left as an Unharvested Control. The pre-harvest period (2 years) was 2010–11 and post-harvest (2 years) was 2012–13. Water table elevation at the upslope and downslope position, subsurface runoff, and soil moisture were measured between May–November. Mixed effect statistical models were used to compare both the before-after and “control” treatment ratios (ratios between harvested hillslopes and the Unharvested Control hillslope). Subsurface runoff significantly increased (p < .05) at both harvested hillslopes but to a greater degree on the Biomass Left hillslope. Greater subsurface runoff volumes at both harvested hillslopes were driven by substantial increases during fall, with additional significant increases during summer on the Biomass Left hillslope. The hydrological connectivity, inferred from event runoff ratios, increased due to harvesting at both hillslopes but only significantly on the Biomass Left hillslope. The winter harvest minimized soil disturbance, resulting in no change to the effective hydraulic conductivity distribution with depth. Thus, the observed hydrological changes were driven by increased effective precipitation and decreased evapotranspiration, increasing the duration that both harvested hillslopes were hydrologically active. The harvesting of residual biomass appears to lessen hydrological connectivity relative to leaving residual biomass on the hillslope, potentially decreasing downstream hydrological impacts of similar forestry operations.  相似文献   

Colorful Planets,Cometary Tail,and Nuclear Winter     

A.?S.?GinzburgEmail author 《Izvestiya Atmospheric and Oceanic Physics》2018,54(3):265-274

After the 1967 flight of the Venera-4 automatic interplanetary station (AIS), the study of planetary atmospheres became one of the key areas of scientific research at the Institute of Atmospheric Physics (IAPh), USSR Academy of Sciences. The goal of the Venera-4 mission was to deliver a descent vehicle into the atmosphere of Venus and study the physical parameters and chemical composition of the atmosphere. IAPh scientists published several articles on the analysis of these measurements, with A.M. Obukhov himself directly involved in writing some of these articles. An analysis of the properties of the Martian surface and atmosphere from the data of the Mars-2, Mars-3, and Mariner-9 missions in November–December 1971 set the stage for a series of studies on the atmospheric effects of Martian dust storms. Based on the study of the atmospheres of terrestrial planets, IAPh scientists developed a simple model for the nuclear winter phenomenon, i.e., a long-term cooling due to massive fires caused by nuclear explosions. Obukhov took a keen interest in this subject and participated in the publication of review articles on the possible atmospheric and climatic effects of a nuclear war. In another publication, Obukhov and his coauthors provided a theoretical analysis of the possible causes for the tail from Halley’s comet ripping off, as was observed in January 1986. The present article gives a brief overview of the IAPh works on Solar System research and on the possible consequences of a nuclear conflict, which were published in the 1960s–1980s while Obukhov was alive.  相似文献   

Projections of the impacts of climate change on the water deficit and on the precipitation erosive indexes in Mantaro River Basin,Peru     

Sly C. Wongchuig  Carlos R. Mello  Sin C. Chou 《山地科学学报》2018,15(2):264-279

Projections of climate change are essential to guide sustainable development plans in the tropical Andean countries such as Peru. This study assessed the projections of precipitation and potential evaporation, rain erosive potential, and precipitation concentration in the Mantaro River Basin, in the Peruvian Andes, which is important for agriculture and energy production in Peru. We assumed the Intergovernmental Panel on Climate Change (IPCC) A1B greenhouse gas emission scenario and simulated the global climate change by the HadCM3 global climate model. Due to the steepness of the mountain slopes and the narrowness of the river valley, this study uses the downscaling of the global model simulations by the regional Eta model down to 20-km resolution. The downscaling projections show decrease in the monthly precipitation with respect to the baseline period, especially during the rainy season, between February and April, until the end of the 21st century. Meanwhile, a progressive increase in the monthly evaporation from the baseline period is projected. The Modified Fournier Index (MFI) shows a statistically significant downward trend in the Mantaro River Basin, which suggests a possible reduction in the rain erosive potential. The Precipitation Concentration Index (PCI) shows a statistically significant increasing trend, which indicates increasingly more irregular temporal distribution of precipitation towards the end of the century. The results of this study allow us to conclude that there should be a gradual increase in water deficit and precipitation concentration. Both changes can be negative for agriculture, power generation, and water supply in the Mantaro River Basin in Peru.  相似文献   

Steady Lateral Growth of Three-Dimensional Particle Laden Density Currents     

M.R. Moossavi Hekmati  M. Najafi  N. Ashrafi khorasani 《中国海洋工程》2018,32(4):467-475

In this paper the steady lateral growth of three-dimensional turbulent inclined turbidity current is investigated. To simulate the current, an experimental setup is developed to analyze the turbidity current for different regimes in the particle laden density currents environment. The Buckingham's π theorem together with a dimensional analysis is implemented to derive the appropriate non-dimensional variables. The experimental results were normalized and plotted in the form of non-dimensional graphs from which a theoretical model is developed and analyzed. Based on the results obtained for the steady lateral growth, three different regimes, namely, inertia-viscous one as the first regime, buoyancy-viscous and gravity-viscous as the second and third regimes are distinguished within the current.In these regimes, the force balance is between the driving and resisting forces. Namely, in the first regime, the force balance is between the inertia and viscous forces, in the second regime, the buoyancy and viscous forces, and in the third regime, gravity and viscous forces are balanced. The experimental results indicate that the lateral growth rate in the first regime is smaller than that in the second and third regimes due to the magnitude and type of the forces involved in those regimes. According to the graphical results, the three different lateral growth rates appear when the normalized current length is smaller than about 3, between about 3 and 10, and larger than about 10. In those regions,the slopes of the data are different with respect to one another.  相似文献   

Projectile preservation during oblique hypervelocity impacts          下载免费PDF全文

Ronald T. Daly  Peter H. Schultz 《Meteoritics & planetary science》2018,53(7):1364-1390

Impact angle plays a significant role in determining the fate of the projectile. In this study, we use a suite of hypervelocity impact experiments to reveal how impact angle affects the preservation, distribution, and physical state of projectile residues in impact craters. Diverse types of projectiles, including amorphous silicates, crystalline silicates, and aluminum, in two sizes (6.35 and 12.7 mm), were launched into blocks of copper or 6061 aluminum at speeds between 1.9 and 5.7 km s⁻¹. Crater interiors preserve projectile residues in all cases, including conditions relevant to the asteroid belt. These residues consist of projectile fragments or projectile-rich glasses, depending on impact conditions. During oblique impacts at 30° and 45°, the uprange crater wall preserves crystalline fragments of the projectile. The fragments of water-rich projectiles such as antigorite remain hydrated. Several factors contribute to enhanced preservation on the uprange wall, including a weaker shock uprange, uprange acceleration as the shock reflects off the back of the projectile, and rapid quenching of melts along the projectile–target interface. These findings have two broader implications. First, the results suggest a new collection strategy for flyby sample return missions. Second, these results predict that the M-type asteroid Psyche should bear exogenic, impactor-derived debris.  相似文献   

Frontiers and route‐ways from Europe: the Early Middle Palaeolithic of Britain          下载免费PDF全文

Nick Ashton  Claire R. E. Harris  Simon G. Lewis 《第四纪科学杂志》2018,33(2):194-211

  相似文献   

Solution for a plane strain rough‐walled hydraulic fracture driven by turbulent fluid through impermeable rock          下载免费PDF全文

N. Zolfaghari  E. Dontsov  A. P. Bunger 《国际地质力学数值与分析法杂志》2018,42(4):587-617

The impact of turbulent flow on plane strain fluid‐driven crack propagation is an important but still poorly understood consideration in hydraulic fracture modeling. The changes that hydraulic fracturing has experienced over the past decade, especially in the area of fracturing fluids, have played a major role in the transition of the typical fluid regime from laminar to turbulent flow. Motivated by the increasing preponderance of high‐rate, water‐driven hydraulic fractures with high Reynolds number, we present a semianalytical solution for the propagation of a plane strain hydraulic fracture driven by a turbulent fluid in an impermeable formation. The formulation uses a power law relationship between the Darcy‐Weisbach friction factor and the scale of the fracture roughness, where one specific manifestation of this generalized friction factor is the classical Gauckler‐Manning‐Strickler approximation for turbulent flow in a rough‐walled channel. Conservation of mass, elasticity, and crack propagation are also solved simultaneously. We obtain a semianalytical solution using an orthogonal polynomial series. An approximate closed‐form solution is enabled by a choice of orthogonal polynomials embedding the near‐tip asymptotic behavior and thus giving very rapid convergence; a precise solution is obtained with 2 terms of the series. By comparison with numerical simulations, we show that the transition region between the laminar and turbulent regimes can be relatively small so that full solutions can often be well approximated by either a fully laminar or fully turbulent solution.  相似文献   

Correction to: Structure of the Transition Region and the Low Corona from TRACE and SDO Observations Near the Limb     

Alissandrakis  C. E.  Valentino  A. 《Solar physics》2019,294(10):1-12

  相似文献